Solitary waves of Bose-Einstein condensed atoms confined in finite rings 
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Motivated by recent progress in trapping Bose-Einstein condensed atoms in toroidal potentials, we 
examine solitary-wave solutions of the nonlinear Schrodinger equation subject to periodic boundary 
conditions. When the circumference of the ring is much larger than the size of the wave, the density 
profile is well approximated by that of an infinite ring, however the density and the velocity of 
propagation cannot vanish simultaneously. When the size of the ring becomes comparable to the 
size of the wave, the density variation becomes sinusoidal and the velocity of propagation saturates 
to a constant value. 

PACS numbers: 05.30.Jp, 67.85. Hj, 67.85. Jk 
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I. INTRODUCTION 

The problem of solitary-wave solutions of the one- 
dimensional nonlinear Schrodinger equation has been 
studied extensively, starting with the work of Zakharov 
and Shabat [l| and of Tsuzuki 2]. In addition, Lieb 
Q has examined the problem of a Bose gas in strictly 
one dimension with periodic boundary conditions and a 
repulsive interaction, and he derived the excitation spec- 
trum. As demonstrated in that study, the spectrum con- 
sists of the usual Bogoliubov branch as well as a second 
branch. This branch was later identified as corresponding 
to solitary- wave excitation 0, @ • 

Bose-Einstein condensed atoms interacting via a con- 
tact potential are described by a nonlinear Schrodinger 
equation and thus offer a very suitable system for the 
study of solitary waves @. Solitary waves have been 
created in quasi-one-dimensional traps with a very weak 
harmonic potential in the long direction |3] and in pan- 
cake traps (8|. Recently, it also became possible to con- 
fine atoms in toroidal trap s (9l-[l3j in which persistent 
currents were observed |13j . 

Based on this experimental progress, it is natural to 
consider the possible creation of solitary waves in such 
toroidal traps. This is the subject of the present study. 
In the following, we examine the properties of solitary 
waves, including their density, propagation velocity, and 
phase (and thus the superfluid velocity), with periodic 
boundary conditions. We note that Carr et al. [ij] have 
previously considered the special case of static periodic 
solitary-wave solutions for a fixed particle number and a 
length. 

For simplicity we consider in the present study a very 
tight toroidal trap, neglecting the transverse degrees of 
freedom, and thus consider the purely one-dimensional 
motion of atoms on a ring. The main question to be con- 
sidered here is the behavior of solitonic solutions as the 
length of the period L is varied for fixed N /L, where N 
is the number of particles. The relevant dimensionless 
parameter is the ratio of L and the coherence length, ^, 
which we choose as that corresponding to the maximum 



density. According to the value of L/^ we distinguish be- 
tween "large" rings, where L/^ 3> 1, and "small" rings, 
where < 1. In Lieb's study [||, the limit iV — > oo 
and L — oo with N/L finite was considered. This cor- 
responds to the limit of large rings in the present study. 
We note that the ratio {S./L)'^ also describes the ratio 
between the typical zero-point energy, ?i^/(2ML^), and 
the typical interaction energy, /(2M^'^). Finally, we 
restrict ourselves to the case of repulsive effective inter- 
atomic interactions. 

As we will show below, in the limit of large rings, the 
wave is localized within a length scale determined by the 
coherence length, and the density is constant over the 
rest of the system, in agreement with known results for 
an infinite system. On the other hand, depending on 
the phase difference at the edges of the system (which 
is either zero or some integer multiple of 27r), two differ- 
ent situations can arise. For zero phase difference, the 
minimum possible value of the density is zero, and the 
velocity of propagation varies between a lowest nonzero 
value, which scales as ^/L and essentially corresponds 
to the superfluid "drift" velocity, and a maximum veloc- 
ity which is set by the speed of sound. In the case of a 
nonzero phase difference, the velocity of propagation can 
be zero, but the corresponding density cannot vanish, 
scaling like (^/L)^. Clearly, when N oo and L ~> oo 
with N/L finite, we find agreement with known results. 

In the opposite limit of small rings, the density be- 
comes sinusoidal with a length scale given by L itself. 
Furthermore, the velocity of propagation saturates to a 
constant value, which scales like 1/L, and the minimum 
density can become zero. 

In what follows we first formulate our problem in 
Sec. II. Then, in Sec. Ill we give the general solution. Sec- 
tion IV examines certain limiting cases, starting from the 
linearization of the problem and then turning to the lim- 
its of large and small rings. In Sec. V we present a general 
result for the minimum value of the velocity of propaga- 
tion of the wave, and finally we summarize our results 
and present our conclusions in Sec. VI. 
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II. FORMULATION OF THE PROBLEM 

Assume that we have a ring of circumference L which 
contains N atoms of mass M . Within the mean- field 
approximation, the order parameter $(a;,i) satisfies the 
equation (in the range —L/2 < x < L/2), 



2M dx^ 



Uo\^{x,t)\''^{x,t), (1) 



with j!^l^/2 \ ^\'^ dx — N and Uq > 0, because of the repul- 
sive interatomic interactions that we consider. Assuming 
that ^{x,t) — 4'(z)e~*''*/'', where fi is the chemical po- 
tential and z = X ~ ut, 



oz 



(C/o|*|'-m)*- 



(2) 



Here u is the velocity of propagation of the assumed 
traveling- wave solution. Setting ^^iz) = y/n{z)e^'^'^^\ the 
superfluid velocity v = {h/M)ip'{z) is given by the conti- 
nuity equation 



v{z) 



]Vr 2APn{z)' 
There is also an Euler equation for the density, 

2^ LVV J \ T -r 3 2 2 



(3) 



(4) 



where Ci and C2 are constants of integration. 

The desired solution has a vanishing derivative at the 



points = and ±L/2 where n 



and n = n-n 



respectively. This allows us to eliminate C2 from Eq. (jH) 
and to determine that 



This leads to 



min"'max 



1. 



.(5) 



— (n')2 = 2J7o(nmax - n){n - + no), (6) 

with no = (2/z -|- Mu^)/Uo — n„-^i^^ — n^^^^- In the case of 
an infinite system, /i = n-maxC^ f^nd n^n^Uo = Mv?, in 
which case Eq. ([6]) then reduces to the correct and known 
form of the differential equation. 



III. THE GENERAL SOLUTION 

The general solution of Eq. ([6|) which is periodic in the 
interval between —L/2 and i/2 is [ij] 



n{z) = riniin -\ 

m 



^min) 



l-dn^ /^2^ 



_ . 2i'^K{m)z^ 

^min ~l" (^max ^minj Sll 



(7) 



where dn(a;|m) and sn(a;|m) are Jacobi elliptic functions 
and K{m) is the elliptic integral of the first kind. Since 
sn(a; — 0|to) — and sn(a; = ■izK{m)\m) — 1, we see that 
n(z = 0) = riniin and n{z — ±L/2) — rimax- In order for 
Eq. ([7]) to be the solution of the differential equation of 
Eq.®, 



K{m) 



2V2 ^ 



1/2 



(8) 



where ^ is the coherence length corresponding to a den- 
sity n-niax, /{2M(^) = nmaxC^O, and A = nnlin/?^max■ 

Similarly, /i and Ci are given as 



A* = 



2 









( 2 - — ) n„iin + 




1 '^max 


\ m J 




m J 



and 
SAP 



Uo 



''mm "'max 



1-1 



(9) 



(10) 



Since K{m) — >-7r/2asm— >0 and K{m) — >■ 00 as m — >■ 1, 
m ranges between and 1. Clearly the solution of Eq. ([7]) 
can be immediately extended to describe k, identical soli- 
tary waves on a ring of circumference kL. 

The velocity of propagation, u, can then be determined 
by integrating the continuity equation, Eq. over one 
period. 



L/2 



V dz 



L/2 



-[0(L/2)-(/)(-L/2)] 



uL 



n^Ci /"^/^ dz 

2AP J-L/2^' 



(11) 



Assuming that 0(L/2) — 0(— L/2) = 2TTq with q an integer 
as a consequence of the periodic boundary conditions, we 
solve the above equation in terms of Ci and combine it 
with Eq. ([S]) to obtain 



u _ 2^/2'Kqi 1 
c ^ L L 



A + 



1-A 



-1 1/2 „L/2 



L/2 



\/ ^min^max 

n(2;) 



dz. 



(12) 



Here c is the speed of sound of a homogeneous gas of 



density 



JJq. We remark from Eq. (fTTj) 



that C\ and u have opposite signs when g = 0. 

While the parameterization of n(z) that we adopt 
above is natural mathematically, it is of greater physi- 
cal interest to explore solutions for fixed iV, L, and [/o 
as a function of, e.g., nmax- To this end, it is convenient 
to pick a value of rimax and m, and determine Jimin from 
Eq. ([8]). It is clear from its definition that nmin must be 
greater than or equal to zero for physically meaningful 
solutions. The value of m can then be adjusted to satisfy 
the normalization condition on ni^z). 
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Equations ([7]), ([8]), and (|T2)) are consistent with Bloch's 
theorem [l^. Because of the periodic boundary condi- 
tions imposed, the order parameter 'i^qiz) of the soli- 
tary wave corresponding to the branch with quantum 
number q is given by that obtained for q — hy excit- 
ing the center of mass motion, ^q(z) ~ e^^''^/-^\I>o(z). 
These two solutions have the same density, but there 
is a difference in the velocity of propagation so that 
[u{q) — u{q = 0)]/c — 2\/2nq^/ L, in agreement with 
Eq. (|12l) . Similarly, the energy spectrum consists of a pe- 
riodic part plus an envelope function which results from 
the energy of the center of mass motion. 



IV. LIMITING CASES 
A. Linearization of the problem 

When rimin — ^ ?^inax, one can obtain equations for the 
force and for continuity by linearizing the problem in 
small deviations of the density and the phase from the 
homogeneous solution — -Roe"*", 

q>^{Ro- 5R) e*('^o+<50) ^ ^ {l-a + idff)), (13) 

where a = SR/Rq and 0o ~ 0. When the two are com- 
bined, we obtain 



—a'" + 2{Mu' - n^,^Uo)a' = 0. 



(14) 



Thus, a' — C sin(fcz -I- 4>q) with C and 0o constants. 
Periodicity requires that kL = 2^1 with I an integer. 
Here, we set ^ = 1 so that 



(7) 



2 £2 



c 



In the limit of "large" rings, i.e., L ^ ^, 



u n^e 



In the opposite limit of "small" rings, L <C ^, then 
c L \ 47r2^2 



(15) 



(16) 



(17) 



We note that Eq. ([T5|) can also be derived directly from 
Eq. (|12p provided that one considers a fixed length L and 
then takes n-a^m T^-max, which implies that m — > in 
Eqs. and ©. 



B. Limit of large rings 

We now consider the limiting forms of the solution 
appropriate for large and small rings. For large rings, 
L ^ ^, and the general solution of Eq. ([7]) reduces to 

n{z) = Umin + (j^max - nmin)tanh^ , (18) 



where C = ■\/2C/(l — "min/'^max)"'^^^- This is in agreement 
with the known result, plus corrections of order e~^/'». 
In this limit, m ^ 1 with 



1 — 16 exp 



V2^ 



(1-A) 



1/2 



(19) 



and also Kirn) « ln(4/ ^/l — m). From Eq. (fT2l) we find 
that 



u 

- w 2V27r^ 
c L 



L 



tan 



1 tanh 



(20) 



which reduces to the familiar result u/c — ^ ?^min/?^max 
in the limit of an infinite ring. 

Figure 1 shows u/c as a function of n„iax for the pa- 
rameter choice N jL — 1, L — 20, Uq = 0.6, and q = 
(with = M 1) so that L/^ > 20. Exact results 
(dots) derived from Eq. pT|) . and values obtained using 
the approximation of Eq. pop are given. The agreement 
is remarkably good except in the limit Umax '^min- In 
this limit, Eq. ([20| yields u/c= 1, which does not contain 
the higher-order quadratic term seen in Eq. (|16p . In the 
opposite limit, where rinun ~> 0, Eq. ((20|) gives the exact 
result, as shown in Sec. V. 

In the limit of small rinun Eq. ([^0]) reduces to 



^«^(l + 2,) + ./^, 

C L V ?^max 



(21) 



Consider first the choice q = 0. In this case nmin can 
be arbitrarily small in analogy with a "dark" solitary 
wave in the L oo limit. The velocity of propagation is 
bounded from below by Umin/c — V2£,tt/L, which 
corresponds to the minimum value oiu/c shown in Fig. 1. 
Equivalently, 



h TT 

ML 



(22) 



except in the immediate vicinity of the solitary wave. 
This result is easily understood. As Eqs. ^ and ([TU]) 
imply, when nmin = 0, Ci = and the superfluid velocity 
is independent of position and equal to the velocity of 
propagation u. In the limit of a dark solitary wave, the 
phase 0(z) varies linearly, except in the limited region of 
small density, where the phase develops a discontinuity of 
— TT. In order to satisfy the condition (p{L/2) = (j){—L/2), 
(j){z) must have the small slope w tt/L over a large length 
w L. Thus, as Eq. (1^^ indicates, u cannot be zero. 
In the case q ^ 0, u can vanish, but when u vanishes, 



27T^e 



(l + 2g)2 



(23) 



This result is also easy to understand. Let us set 5 = — 1 . 
In this case Eq. implies that v{z) oc l/n{z). Away 



4 



1.2 




1 1.04 1.08 1.12 1.16 

"max/"o 

FIG. 1: The velocity of propagation ujc in "large" rings, as a 
function of nmax/no, where no = N/ L, with no = 1, L = 20, 
Uq = 0.6, and q — 0, with ft = M = 1. Here the dots give the 
exact result, Eq. and the solid line gives the approximate 
result of Eq. ([20}. In this figure L/^ > 20. 



10 




4 



2 ■ 

I ■ ■ ■ ■ — 
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FIG. 2: The velocity of propagation u/c in "small" rings, as a 
function of nmax/no, where no = N/L, with no = 1, L = 0.5, 
Uo = 0.6, and q ~ 0, with h — M = 1. Here the dots give the 
exact result, Eq. and the solid line gives the approximate 
result of Eq. In this figure L/^ < 0.775. 



from the center of the wave, v{z) is constant and equal to 
{'h/M)(Tr / L) since the density is constant, but the phase 
(l){z) is equal to it{z/L). The total accumulation of the 
phase over L is thus approximately tt, while another tt is 
needed in order to create the total phase difference of 2n. 
This is provided from the region around the wave where 
v{z) K, {'h/M){n I L) ny^s.^/n{z), within a length scale of 
order 'T-min/'^max- This is the length scale of the dom- 
inant integration interval of v{z) around its singularity. 
Since (f>' is ~ (7i'/i)(n,nax/"-min) in this regime, and in 
order for the integral of (l)'{z) to be equal to tt when we 
integrate it around a length scale of order ^ \/ 'T.min/'^max, 
we find that 



4>'{z) dz 



^ ^max 



(24) 



This implies that nmin/»^max ~ (C/^)^- 

These results agree with those of Carr et al. (T^]. In 
that study only static solitary waves were considered, and 
(as mentioned in Ref. this is not possible for zero 
phase difference. Their choice of ^/L = 1/25 leads to 
'T-min/'^max ~ 0.032 according to Eq. (j23|) . which is in 
good agreement with the numerical results of Fig. 4 of 
this reference. 

It should also be noted that Eq. ([23| is valid for a single 
solitary-wave solution. For a two-solitary wave solution 
it is possible for them to be both static and dark (i.e., to 
have a node in the density). 

The maximum possible value of u is achieved in the 



limit of sound waves n-n 



where 



l + 2N/27rg| + ^, 



(25) 



in agreement with Eq. p^ . Finally, we note that 
the normalization condition takes the form N/L = 

1 - 2v^(C/L) (1 A)^''^] in this limit. 



C. Limit of small rings 

For small rings, L ^ f , the general solution of Eq. ([7]) 
reduces to 



n{z) 



I sm 



(26) 



plus corrections of order m (i.e., of order (i/^)^) since 
TO sa (1 - A)(LV27r2^2) ^(^^ ^ (7r/2)(l + m/4) in 
the limit to — !■ 0. Turning to the velocity of propagation 
of the wave, Eq. p2)) implies that in the limit of small 
rings 



u 
c 



L 



(1 + 2?) 



L 



2\/27r^ 



(27) 



One can identify the term proportional to 1 -I- 2(7 as the 
difference (£^,+1 - Eq)/5p, where Eq = h^q'^ /{2MR'^) is 
the single-particle kinetic energy of a particle in a ring of 
radius R with angular momentum hq, and 6p = h/R. 
This expression coincides with Eq. P7)) . found by lin- 
earization, when q = and rimin — nmax- Here, rimm 
is not bounded from below as it was in the limit of large 
rings. Remarkably, u saturates to a constant value in- 
dependent of ?T.i„in/nmax whcu L / 0. Tliis happens 
because 



u 
c 



2V2- 



TTq 



^ 1 fi-xv^^ r'^ 



L 



L 



(1 + 2<7) 



7r/2 A -f (1 - A) sin^ w 
(28) 



in this limit, since the integral is independent of A. In 
other words, in the limit of small rings, the sinusoidal de- 
pendence of the density always gives rise to a phase dif- 
ference of TT for any value of the ratio TT-min/'^max- Since 
the periodic boundary conditions require that the total 
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phase difference has to be 27rg, the additional phase ac- 
cumulation must be (1 + 2g)7r. This fact gives a slope to 
the phase that is w {l + 2q)TT/ L (plus corrections of order 
(i/^)^, and therefore w « (1 + 2q){fi/M){-K / L), in agree- 



ment with Eq. (1281) . Finally, the normalization condition 
is N/ L — (nniin + ?T-max)/2 in this limit. 

Figure 2 shows u/c as a function of rimax for the pa- 
rameter choice N/L = 1, L = 0.5, Uq — 0.6, and q = 
(with h = M = 1), so that L/C < 0.775. Exact results 
(dots) and values obtained using the approximation of 
Eq. (|27l) are given. The agreement is remarkably good 
in both limits rtmax '^mim and rimin — > 0, as Eq. (P7)) 
gives the exact result in both cases, as Eqs. (|17l) and ((29)) 
indicate. 



V. A GENERAL RESULT 

It is interesting to note that, in the limit of a dark 
solitary wave (nmin 0), the velocity of propagation u 
of the wave is given by the same formula. 



(l + 2g). 



or 



b' = -{l + 2q), 



(29) 



(30) 



in the limits of both small and large rings, as indicated 
by Eqs. (I21|) and dST]). In fact, Eq. ^ gives the mini- 
mum possible value of u for any L: A change of variable 
tranforms Eq. (fT^ into 



u _ 2^/27^(^C 
c 

i/Vx 



1 



L 



2K{m) 
dw 



A + 



1 - A 



-,1/2 



m 



-i/Va [1 + (1 - A)w2]Vl - AwVl - Amw2 



(31) 



The integral can be evaluated in the limit A — ?► by 
closing the integration contour with a semicircle in the 
upper half complex plane. The only contribution to the 
integral comes from the pole dX w — i. In this way, one 
obtains Eq. (P^ for any finite value of L. 

In the limit where n^in — > 0, the integral in Eq. (fT^ 
is dominated by the behavior of the density near z — Q 
since sv?{x\m) for a; — )■ 0. As a result, the integral is 
approximately proportional to dz/{X + z^) = tt/a/A. 
As a result, there is a contribution to (j)' which is equal to 
tt/L. For a nonzero q, there is an additional contribution 
of 2TTq/L to (j)', from which Eqs. ((291) and dSO]) follow. 



tures if one imposes the periodic boundary conditions ap- 
propriate for the description of Bose-Einstein condensed 
atoms confined in tight toroidal traps. 

In the limit of large rings, where the size of the wave 
is much smaller than the length of the ring, the density 
is exponentially localized on a length scale determined 
by the coherence length and is constant over the remain- 
der of the ring. Depending on the phase difference at 
the edges of the period, either the minimum density can 
vanish (with the velocity of propagation being nonzero) 
or the velocity of propagation can vanish (with the min- 
imum density being nonzero). However, both cannot be 
zero simultaneously. 

In the limit of small rings, the density is sinusoidal 
within a length scale given by the circumference of the 
ring itself. The velocity of propagation saturates to a 
constant value, which is independent of nmin. For this 
reason, the size of the ring must exceed some minimum 
size in order to create a static wave with u = [14t] . 
The minimum size can be determined from Eq. (j3ip for 
q — —1 and in the limit to — > and A — 1 with to/(1 — 
A) = 3. Under these conditions, u = and L/£^ ~ \/6 7r, 
in agreement with Ref. [14 1. 



In the limit of an infinite ring, the velocity of propaga- 
tion of the solitary wave depends on the ratio rimin/'^max- 
In the opposite limit of small rings, the range of possible 
values of the velocity of propagation becomes narrower as 
L decreases. In the limit of very small L, u saturates to a 
constant value and becomes independent of nmin/?T-max- 

The general result that the minimum value of u is in- 
versely proportional to L may be relatively easy to inves- 
tigate experimentally. 

Finally, we note that as opposed to the case of an in- 
finite ring, which contains only the dimensionless quan- 
tity n„iin/'T-max, & sccoud dimcusionless quantity, L/C, 
appears in the present problem. In particular, in the 
limit of sound waves (where 1 — Jimax/'T-max — >■ 0) and 
large rings (where L/C — >■ oo), the product of these two 
terms appears, and the two limits "compete" . Different 
answers will result depending on the way in which these 
limits are taken. These different answers are not only 
of theoretical interest; they represent different physical 
situations determined by the way that an experiment is 
performed. 
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